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The first characterization of mixed-state entanglement was acliieved for two-qubit states in 
Werner's seminal worlc [Phys. Rev. A 40, 4277 (1989)]. A physically important extension of this 
result concerns mixtures of a pure entangled state (such as the Greenberger-Horne-Zeilinger [GHZ] 
state) and the completely unpolarized state. These mixed states serve as benchmark for the ro- 
bustness of entanglement. They share the same symmetries as the GHZ state. We call such states 
GHZ-symmetric. Despite significant progress their multipartite entanglement properties have re- 
mained an open problem. Here we give a complete description of the entanglement in the family of 
three-qubit GHZ-symmetric states and, in particular, of the three-qubit generalized Werner states. 
Our method relies on the appropriate parameterization of the states and on the invariance of entan- 
glement properties under general local operations. An immediate application of our results is the 
definition of a symmetrization witness for the entanglement class of arbitrary three-qubit states. 

PACS numbers: 03.67.-a, 03.67.Mn 
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Introduction. - Entanglement is the essential re- 
source for many tasks in quantum information process- 
ing [2|, [sj . Therefore, it is desirable to precisely charac- 
terize the entanglement contained in a quantum state. 
While our understanding of pure-state entanglement has 
significantly improved in recent years, entanglement in 
mixed states has remained a notoriously difficult subject, 
despite numerous impressive results, e.g., |4l4l2|. Due to 
the tremendous experimental progress in produc ing and 
controlling multi-qubit entanglement (e.g., Refs. J13l4l7|) 
this has become also a practical problem, as an accurate 
assessment of the experimental results is required. The 



universal tool here are entanglement witnesses [8|, ll8l42C 
Despite its flexibility and success in detecting entangle- 
ment and distinguishing entanglement classes, character- 
ization by means of viritnesses is not always satisfactory. 
Enhancing the quality of entanglement witnesses requires 
improvement in the underlying entanglement theory. 

An important and in general unsolved question of both 
practical and theoretical interest is how much noise ad- 
mixture pure-state entanglement can sustain. Mathe- 
matically, this question can be cast as follows. We con- 
sider a pure state Ji/'me) of N qubits that contains a max- 
imum amount of a certain entanglement type. This state 
gets mixed with the operator ^^2^^ which describes the 
maximally mixed state of TV qubits, serving as a model 
of unpolarized noise: 

1 

2^ 



PWS{P)=P |'0Me)(V'Me| + (1-P)7rw^2« ■ (1) 



The question then is how small p (0 < p < 1) can be 
chosen such that pws (p) still contains a finite amount of 
the considered entanglement. 

For two qubits (A^ = 2), one substitutes IV^me) with the 
Bell state |*^) = -^ (|01) - |10)). Then, pws(p) repre- 



sents the so-called Werner states [l|. Although Werner 
defined them through the symmetry under local unitaries 
U ®U , the generalizations of Eq. ^ to three and more 
qubits are often termed generalized Werner states [J,|2l|. 
Throughout this article, we shall consider three qubits 
{N = 3) and the maximally entangled GHZ state 



|GHZ) = ^(1000) + |111» 



I GHZ, 



(2) 



We shall give a complete characterization of the entan- 
glement in 



Pws(p) = p I GHZ) (GHZ I 



l-p 



ils 



(3) 



and the entire family of states with the same symme- 
try, the GHZ symmetry (see below). After reviewing the 
known results we present a parameterization for this fam- 
ily that allows to deduce the entanglement type for any 
given element of the family. Finally we show that our 
findings can be used as a witness to detect the entangle- 
ment type of arbitrary three-qubit states. 

For two qubits, Eq. ^ gives the standard Werner state 
after replacing It/^me) with the Bell state JVP"). In the 
two-qubit case there is only one type of entanglement, 
and therefore the problem reduces to finding the maximal 
value of p such that the state is still not entangled. It 
can easily be found by computing the concurrence [22|. 

The three-qubit case, however, is more complex. A 
state can either be completely separable, biscparablc or 
tripartite entangled. Moreover, there are two inequiva- 
lent classes of tripartite entanglement, the GHZ type and 
the W type [23|. In the space of density matrices, there 
is a hierarchy of entangled states [5| : the convex hull of 
the VF-type states includes the true W states, the bisep- 



arable and the separable ones while the set of GHZ-typc 
states contains all other classes. 

The thrcc-qubit generalized Werner states pws (p) are 
known to be fully separable if and only if p < Psen_ = 
i Q and biseparable if and only if p < Pbisop = f [l3l ■ 
One aim of this article is to find the value pw such that 
Pws {p) is of W type for p < pw and of GHZ type for p > 
Pw- It turns out that it is advantageous to extend the 
problem to all mixed states which can be written as afhne 
combinations of |GHZ+), |GHZ_) = ^jdOOO) - |111)), 
and the maximally mixed state. 

Parameterization of GHZ- symmetric states. - We 
solve the problem by exploiting its symmetry. The GHZ 
state, and thus also pws{p), is invariant under the fol- 
lowing transformations (and combinations thereof): (i) 
qubit permutations, (ii) simultaneous three-qubit flips 
(i.e., application of cr^ (g) ctj: (X) (T^^), (Hi) qubit rotations 
about the z axis of the form 



U{^U^2) 



-i{4>i+<p2)o-z 



(4) 



Here, a^ and Cz are Pauli operators. We refer to the 
invariance under the operations (i) — (Hi) as GHZ sym- 
m,etry. Except the qubit permutations all those opera- 
tions are local, therefore (and since qubit permutations 
always convert GHZ states into GHZ states) GHZ sym- 
metry operations will never turn GHZ-type entanglement 
into ly-type entanglement or vice versa. 

An important aspect of this symmetry is that for any 
decomposition of pws{p) into pure states there is a GHZ- 
symmetric decomposition of the same entanglement type. 
It is generated by replacing each pure state in the decom- 
position with the equal mixture of all states obtained 
from that former state by applying the symmetry opera- 
tions. 

In order to identify the set of GHZ-symmetric density 
matrices we check the action of the symmetry operations 
on its elements p^. First consider the z rotations (Hi). 
The matrix element pfj. ;,„„ is transformed by operations 
according to Eq. @ into exp[i(i — k — 1 + n)(j)i] exp[i(j — 
k — m + ")'/>2]/3fjfc_im„- Since 0i and (j)2 can take arbi- 
trary values the state remains unchanged only if either 
the matrix element is zero, or if both i — k — l + n = and 
j — k — m + n = 0. Therefore the only non-zero matrix 
elements are the diagonal elements, Pooo.iii and Pin qoo- 
Among these elements, permutation invariance forces the 
diagonal elements to depend only on the number of Is in 
the index. Finally, the invariance under collective bit 
flips implies p^oo^ooo = Piii.in ^nd pgoi.ooi = Piio.iio- 



Moreover, we have pggo m = pfu ooo and thus real off- 



diagonal matrix elements due to hermiticity. Given the 
additional constraint trp^ = 1 we find that a state p^ 
is fully specified by two independent real parameters. A 



possible choice is 

^(P') = I [(GHZ+I pS IGHZ+) - (GHZ_ | p^ |GHZ_)](5) 
2/(pS)=-L[(GHZ+|pS|GHZ^ 



(GHZ_|pS|GHZ_)-i 



(6) 



such that the Euclidean metric in the (x, y) plane co- 
incides with the Hilbert-Schmidt metric on the density 
matrices. The completely mixed state is located at the 
origin. The set of states p^ forms a triangle in the {x, y) 
plane (see Fig. 1). The generalized Werner states ([3]) are 
found on the straight line y = ^x connecting the origin 
with the GHZ state. We call it the 'Werner line'. 




FIG. 1. The convex set of GHZ-symmetric density matrices 
p . The upper corners of the triangle are the standard GHZ 
state |GHZ+), and |GHZ_). Note that these are the only pure 
states. Applying a^ to any one of the qubits changes the sign 
of X. Therefore for properties invariant under local unitaries, 
we have a mirror symmetry about the y axis. 
At the centre of the upper horizontal line there is the separable 
state i(|000)(000| + |lll)(lll|). The points for the pure states 

|001), 14-4-4-), \+)\(j)^), \W+ ) indicate the positions of the 

corresponding symmetrized mixed state (here we have used 
the definitions I ±) = ^(|0)±|1)), |0+) = ^(|00>4-|11)), and 

\W+—) = ^(|4- - -) 4- 1- 4- -> 4- 1- - +))). The solid ma- 
genta line ('Werner line') represents the generalized Werner 
states pws(p) with pscp — | and pbiscp ~ f (see text). The 
two lower dashed Unes are first guesses for the boundaries of 
fully separable ('sep') and biseparable ('bisep') states from 
the known values of pscp and pbiscp, respectively. The upper 
dashed fine represents a first guess for the boundary between 

W and GHZ states as \W^ ) is the W state with the largest 

overlap to |GHZ+) (cf. Ref. Ql). The intersection with the 
Werner line occurs at p = 9/13. 

For any normalized pure state \ip) = (^ooo, • ■ • , V'iii)j 
there exists a corresponding symmetrized state 



p^iip)= / dC/C/|V')(V'|C/^ 



(7) 



where the integral is understood to cover the entire GHZ 
symmetry group, i.e., unitaries U{(f)i,(j)2) as in Eq. ((3|) 



and averaging over the discrete symmetries. The coordi- 
nates of the symmetrized state can be inferred from the 
coefficients -0000 and ^m 



1 



XW = 2 (V'oooV'lll + IpOOQlPlu) 



(8) 
(9) 



Entanglement properties of GHZ- symmetric states. 
After finding and suitably parameterizing the set of GHZ- 
symmetric states we want to determine the entanglement 
class of each state (fully separable, biseparable, W , or 
GHZ). The key idea is that all states in an entanglement 
class are equivalent under stochastic local operations and 
classical communication (SLOCC) [23J, |2J]. Mathemati- 



cally, the corresponding (invertible) local operations are 
represented by the elements of the group GL(2, C). That 
is, applying GL(2,C) transformations to any qubit does 
not change the entanglement class of a multi-qubit state. 

The GHZ-symmetric states of each SLOCC class form 
a convex set. We characterize each set by finding its 
boundary starting from the separable states. Our strat- 
egy to identify these boundaries is the following. We fix 
the y coordinate in the interval —1/(4-^/3) < y < a/3/4 
and then consider all pure states |'0) of the SLOCC class 
under consideration whose symmetrized state p {ip) has 
the chosen y according to Eq. ^. States at the bound- 
ary are the ones with maximum (or minimum for x < 0) 
X coordinate according to Eq. ([8]) for a given y, termed 
Xniax- Mirror symmetry implies Xmin = — a^max, therefore 
we may restrict our discussion to a; > 0. If a;niax(y) does 
not have the appropriate curvature the boundary is given 
by the convex hull of a;max(y)- 

We start with an obvious solution that holds for all 
SLOCC classes. If, for fixed j/, the coefficients of the 
pure state \'4>) can be chosen equal IV'oooP = IV'niP = 
\ (a/32/+ j) the maximum x coordinate is given by Xmax = 
iV'ooollV'iiil = \{V?>y + i), i-e., by the lower edge of the 
triangle of GHZ-symmetric states. 

Now consider the separable pure states |i/)''°p). They 
are equivalent (via local unitaries) to the state |000) 



W 



scp\ _ 



J=l 



A] 



B, 



b; -a, 



|000) (10) 



where \Aj\,\ 
efficients we 



B, 



< 1. 



find IV'^"'^ 



000 1 






i^iP)(i 



l^2P)(l-|A3P) 



oooliv^ml subject to the constraint IV'ooo 



For the moduli of the co- 

= \A,A2A3\ and IV^fi?! = 
Maximizing x = 



sop 1 2 



const leads to 



Tsy 



for y > 0. As this 

function gives a concave boundary (cf. Fig. 2) we use the 
convex hull 



pScp 



V3 



1 



(11) 



which is identical to the first guess from the known result 
Pscp = ^ (cf. Fig. 1). 
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FIG. 2. The SLOCC classes of three-qubit GHZ-symmetric 
states p . The dark blue region shows the separable states 

('sep') with the light blue lines x — ±lj i=y) . Green 

areas represent the biseparable states ('bisep'). The W states 
W (yellow) and the GHZ states 'GHZ' (grey) are separated 
by the curve Eq. (|13[l (red line). The Werner line (magenta) 
crosses that curve at pw ~ 0.6955. Some geometrical aspects 
are noteworthy. The curve p3p nearly (within a few per cent) 
describes a circle about the point p (001). The radius has a 
minimum in the vicinity of the Werner line. Further, it is 
intriguing to note that each SLOCC class shares exactly one 
fourth of the lower edge of the triangle. 



For biseparable pure states |?/j^'^°p) it suffices (due to 
the subsequent symmetrization) to consider local equiv- 
alence to the state |O)(8)|0+). That is, we obtain |i/)'^"'''p) 
by normalizing the vector {Gi^ G2^ G3)\0)^\(j)'^) . Here, 



G, 



Aj Bj 



J = 1,2,3 



denotes an arbitrary GL(2,C) transformation. The dis- 
cussion can be restricted to G2 = G3 because for any 
Ujbiscp^ the two-qubit part can be made permutation 
symmetric by a diagonal GL(2,C)®^ operation without 
decreasing the coordinates a;(V'^'^°P), y(V'^'*'°P) of the cor- 
responding symmetrized state. Maximizing x as before 
yields 



V3 



-y + - 
2 '^ 8 



(12) 



for y > —7= . Again this boundary coincides with the one 
inferred from p°^^'^v = | (see Fig. 1). 

The general pure W state |V'^) is found by normaliz- 
ing (Gi®G2®G3X|001)-h|010)-h|100)). In analogy with 
the separable states, maximization of a; = ["iAoooll^ml 



subject to the constraint IV'oool + I '/'nil = const shows 
that the maximum is reached for Gi = G2 = G3. It 
leads to polynomial equations whose solutions are given. 



for convenience, in parameterized form (with < u < 1) 



„w 



-8v^ 



8(4 -u2) 



y 



V34 



(13) 



where y > -y=. The second derivative of x^^^{y) shows 
that the boundary is indeed convex. This completes the 
characterization of SLOCC classes for GHZ-symmetric 
three-qubit states. 

A particularly interesting point is the intersection of 
the curve (fT3|) with the Werner line j/vvs = 2^- "^^^ 
corresponding parameter vw solves the equation 



1 



-"w 



•Jw 



J-W 



« 



such that -pw = 0.6955427 . . . 

Symmetrization witness. - Although these results 
might seem of purely theoretical interest they have a 
surprising application for arbitrary three-qubit states. 
Suppose p is such a state. The twirling operation in 
Eq. ([7|) generates the corresponding symmetrized state 
p^{p}- The SLOCC class of p cannot be lower in the hi- 
erarchy described in the introduction than that of p^ [p) . 
For example, a W state can be projected by the twirling 
operation Eq. ([7]) onto a W state, a biseparable state or 
a fully separable state, but not onto a GHZ state. Hence, 
the GHZ-symmctrized state p^ (p) can be used to witness 
the SLOCC class of the original state p, simply by reading 
off the coordinates of p^{p) in Fig. 2. These coordinates 
x{p) and y{p) are obtained from the matrix elements of 
P- 

x{p) = X (pooo, 111 + Pill, ooo) 

y{p) = ^ (P000.000 + Pill, 111 - ^ 

We will discuss the optimization of this method else- 
where. 

Summarizing, we have determined exactly the en- 
tanglement properties of an entire family of high- 
rank (mostly eight) mixed three-qubit states with the 
same symmetry as the GHZ state. In particular, we 
have solved the problem for the three-qubit generalized 
Werner state which is a reference for multi-qubit mixed- 
state entanglement. A practically relevant application of 
this result is a simple method for detecting the SLOCC 
class of arbitrary three-qubit states. 
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